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Abstract 

Both a and n are well established from E791 data on D — > Stt and 
Ktttt and BES II data on J/^ cj-tt+tt" and K^K~^it^tt^ . These fits 
are accurately consistent with vrvr and Ktt elastic scattering when one 
allows for the Adler zero which arises from Chiral Symmetry Breaking. 
The phase variation with mass is consistent between elastic scatter- 
ing and production data. Also Colangelo et al. show that crossing 
symmetry and dispersion relations for vrvr elastic scattering demand a 
a pole within 2 standard deviations of the pole fitted to production 
data. Oset and collaborators find similar results using unitarised Chi- 
ral Perturbation Theory. Possible interpretations of a, k, /o(980) and 
ao(980) are explored. 

PACS: 13.75CS, 14.40Cs, 14.40.Ev 
Keywords: mesons, resonances 

1 Introduction 

This is an extended version of a report for the Proceedings of the HadronOS 
conference, where space did not allow full discussion. An extensive collection 
of figures of data is given here. Also technical points of discrepancies between 
analyses are discussed in detail. A section involving these technical points is 
denoted by an asterix. 

This work has been assisted greatly by a working group on Scalars, es- 
tablished at HadronOS. I am grateful to the large number of contributors to 
this working group. Opinions have frequently differed, so the responsibility 
is mine for trying to arrive at a consensus; at several points I need to explore 
conflicts of opinion and sometimes try to reach a conclusion. 
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2 The cr pole 



Early evidence for the a pole arose from elastic scattering data. Markushin 

and Lochcr [1] summarise many determinations. Renewed interest was sparked 
off by E791 data on (7r+7r^)7r+ [2]. The tttt mass projection, shown in 

Fig. 1, has a low mass peak which was fitted by a pole shown in the first entry 
of Table 1 below. Their fit assumed a Breit-Wigner resonance with r(s) oc 
where p{s) is Lorentz invariant phase space 2k /^/s — ^1 — Am^/s and 
k is centre of mass momentum. This choice for r(s) will later be shown to 
be inappropriate, but correcting it to a better form introduces only changes 
of detail. Oiler has refitted the data in a way consistent with Chiral Pertur- 
bation Theory; there is some increase in the width of the pole [3]. 




Figure 1: The tttt mass projection of E791 data for —>■ tt~^{tt tt"*"), (a) 
without, (b) with a in the fit. 



Higher statistics data from BES II [4] on J/^ — > a;7r+7r~ are shown in 
Fig. 2. They are dominated by /2(1270), 6i(1235) and a, which is clearly 
visible as a flat band along the upper right-hand edge of the Dalitz plot in 
(b). The 0+ contribution is shown in (e). There is a marginal 2 standard 
deviation contribution from /o(980). In order to test the sensitivity of the 
extrapolation to the pole, four types of parametrisation were tried. All are 
consistent with an average pole position M = (541 ± 39) — i(252 ± 42) MeV. 

3 How to parametrise the a 

Fig. 3(a) shows the intensity of tttt elastic scattering v. mass. Why is there 
no low mass peak like that in production data of Fig. 2? The explanation was 
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Figure 2: BES II data for J/\I' cutt+tt^. (a) The lo mass peak showing 
background from 7r+7r^7r+7r^7r^; (b) Dalitz plot; (c) vrvr mass projection; the 
histogram shows the fit and the hatched area the experimental background; 
(d) 101^ mass projection; the dashed histogram shows the 61(1235)71 contri- 
bution (two combinations); (e) the a contribution (hatched) and the full 0"'" 
contribution including /o(980); (f) 2+ contribution. 
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given in 1965-6 by Adler and Weinberg [5]. They proposed that massless tt 
of zero momentum have zero elastic amphtude. If the / = S-wave nn — > mr 
amphtude is expanded as a power series / = arn^ + bk'^, consistency between 
s, t and u channels requires / oc (s — 0.5m^) and a zero at the Adler point 
= 0.5m^. Fig. 3(b) shows the result of dividing Fig. 3(a) by (s — Syi)^. 
Instantly one sees a resemblance with the a peak of Fig. 2. So the solution 
to the puzzle is that the matrix element for tttt elastic scattering is strongly 
s-dependent: a situation unlike most other resonances. 




Figure 3: (a) The tttt intensity in elastic scattering, (b) with the Adler zero 
divided out. 



Let us write the elastic a amplitude as 

Njs) _ N^ijs) 

D{s) - s - iNtotis)' ^ ' 

Here Nei{s) is real for s > 0. The phase variation comes purely from -D(s); 
this denominator is universal for all processes involving a tttt pair. [That is 
the assumption on which the Particle Data Tables are based]. For elastic 
scattering, the Adler zero in N{s) nearly cancels the a pole for low masses. 
However, the numerator N{s) is not universal; it is quite different for pro- 
duction processes, where the left-hand cut is distant. Later, it will be shown 
that E791 data for {K~TT~^)7r~^ require N{s)prodn — 1 within errors. 

The production amplitude will therefore be written 

TSLn-^/D{s), (2) 
where A is a complex coupling constant. 
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Figure 4: Contours of intensity for (a) production, (b) elastic scattering. 



Fig. 4 sketches contours of constant intensity for (a) production, (b) 
elastic scattering. In the latter case, the effect of the Adler zero is to suppress 
the intensity near s = 0. The elastic phase shift on the real s axis (where 
experiments are done) reaches 90° only at M > 900 MeV, far above the 
position of the pole. It is this feature which confuses many people. However, 
the phase varies rapidly off the real axis because the width increases with s 
and because of consequent effects of analyticity. 

A simultaneous fit is made to BES II data, Kg4^ [6] and Cern- Munich data 
[7] using the empirical form: 



N{s) = M{s - 0.5ml) exp[-(s - M'')/A]{l + Ps)p^^{s) + Mr^^{s). (3) 



The exponential is required by nn elastic data to cut off N{s) above 1 GeV. 
The term Mr4^ accounts for An inelasticity above ~ 1 GeV, but has little 
effect on the a pole. Note from Fig. 2(e) that the a intensity fitted to BES 
data falls to a low value above 1 GeV. The a pole is therefore distinct from 
the broad /o(1535) fitted by Anisovich et al. [8]. It is also distinct from the 
broad pole fitted around 1 GeV by Au, Morgan and Pennington [9]. 

Below the KK threshold, the elastic amphtude must follow the unitarity 
circle. There are small contributions from the low mass tails of /o(1370), 
/o(1500), etc. and a contribution from /o(980) which is large around 1 GeV. 
These are included by writing 




l)/(2i), 



(4) 



where 5* is the S'-matrix riexp{2i6) = 2i(l + T) for each individual resonance. 
Below the KK threshold, = 1 for all amplitudes, so phases due to individual 
resonances add. 

For production, there are hundreds of open channels for D and J/^ de- 
cays. Within individual channels, unitarity plays a negligible role. In the 
standard isobar model, amplitudes are added using a complex coupling con- 
stant A = 5fexp(i0o) ioT each amplitude; the phase (pQ is constant over the 
whole Dahtz plot. For J/^^ — > Loa, 0o is the phase of auj elastic scattering at 
the J/^ mass; it is unknown, so g and 0o need to be fitted to the data. 

The Ke4 data are available up to 380 MeV and there is then a gap in 
elastic scattering data until 560 MeV, where Cern-Munich data begin. The 
a pole lies in the mass range where there are no elastic data. Although this 
gap may be bridged by using dispersion relations, the production data are 
obviously important in filling the gap directly. 

4 The phase of the a 

In the BES data, the bin channel contributes 41% intensity and auj 19%. 
There are strong interferences between them which determine the phase vari- 
ation of the a with mass. The data have been divided into bins 100 MeV 
wide from 400 to 1000 MeV. Fig. 5 shows phases for individual bins, keeping 
magnitudes fixed at values from the global fit: this achieves optimum accu- 
racy, since there is noise in the magnitude in individual bins if it is set free. 
The full curve shows the phase from elastic scattering data. The agreement 
with the bin-by-bin fit shows that the Breit-Wigner form of eqn. (1) indeed 
accounts for the phase in all data. The imphcation is that the data at low 
mass may be described by a single resonance, except for well understood 
contributions from /o(980), /o(1370) and /o(1500). It is not necessary to 
include a background amplitude. Another way of stating this is that if any 
background is present, it is the same in both elastic scattering and produc- 
tion and may be absorbed algebraically into the parametrisation of the a 
amplitude. 

Within the scalar working group, it has been suggested that phases could 
be modified by rescattering of pions from a decay off the uj; although this 
must happen at some level, its contribution seems to be below the level of 
experimental errors at present. 
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Figure 5: The phase of the a amphtude in mass bins 100 MeV wide, compared 
with the global fit (full curve), a Breit-Wigner of constant width (dashed) 
and a Breit-Wigner with V oc Ptttt (dotted). 



The dotted curve of Fig. 5 shows the phase variation if the amphtude for 
the a is expressed in the form F oc Ptttt- This gives a marginally poorer fit, but 
cannot be distinguished cleanly using the BES data. Unfortunately, the hi 
band runs off the corner of the Dalitz plot and does not interfere significantly 
with the a below 400 MeV. 



It is important to realise that there is a large background of theoretical work 
on elastic scattering and related processes like tttt — > KK. 

Colangelo, Gasser and Leutwyler [10] have made a precise determination 
of the a pole from elastic scattering and K^^^ data without using production 
data. They use crossing symmetry to calculate the amplitude on the left-hand 
cut (s < 0) from amplitudes at s > 4m^ (right-hand cut). For a given s and 
t, the two are related by a simple isospin matrix; the tttt S-wave amplitudes 
is then evaluated from an integration over the appropriate range of t (and 
likewise for P and D- waves, including appropriate Legendre polynomials). 
As a further constraint, they use fixed t dispersion relations (Roy equations). 



5 Theory 
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State Reference Data Pole position 



(MeV) 

(7 [2] L)+ (7r+7r-)^+ (489 ± 26) - i(173 ± 26) 

[3] D+ (7r+7r-)7r+ 470 - i220 

[4] J/^ ^ a;(7r+7r-) (541 ± 39) - i(252 ± 42) 

[10] TTTT ^ TTTT (470 ± 30) - i(295 ± 20) 

[13] TTTT TTTT 445 - i221 

[25] TTTT ^ nn (470 ± 50) - i(285 ± 50) 

[27] D+ {n+n-)n+ (455 ± 36) - i(190 ± 36) 

K [30] D+ ^ {K-71+)ti+ (721±61) -i(292± 131) 

[32] J/^ ^ K+7r-K-7r+ (760 ± 41) - i(420 ± 75) 

[31] J/* ^ K+7T-K-n+ (841 ± 82) - i(309 ± 87) 

[36] Kn Kn (722 ± 60) - i(386 ± 50) 

here all 750l-^^ - ^(342 ± 60) 

[3] D+ {K-7r+)Tr+ 710 - i310 

[12] Kn Kn (770 - i(250 - 425) 

[14] Kn Kn (708 - i305) 

[15] Kn Kn (753 - i235) 

[24] Kn^Kn (594 ± 79) - i(362 ± 322) 

/o(980) [37] J/^ ^ (/)7r+7r- (998 ± 4) - i(17 ± 4) 

[12] nn nn and KK 994 - il4 

ao(980) [39] pp rjnn and Lorjn^ (1036 ± 5) - ^(84 ± 9) 



Table 1: Summary of pole positions. 



This extends greatly the range of s in which the S-wave is known and fixes the 
/ = and 2 scattering lengths accurately. In this work, Re is accurately 
determined and contains a very clear Adler zero at s = 0.45m^^ (displaced 
by a tiny amount from Weinberg's prediction because of higher powers of k). 
The pubhshed work fits nn masses up to 800 MeV and demands a a pole at 
M = (470±30) -i(295±20) MeV, i.e. about 2 standard deviations from the 
BES experimental value. Further work is in progress, to extend the fitted 
mass range to 1150 MeV, above /o(980) and the KK threshold. 

In a series of papers [11-16], Oset, Oiler, Pelaez and collaborators have 
fitted data using 'unitarised' Chiral Perturbation Theory. In outline, they 
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use ChPT for lowest order and take rescattering from the next order. They 
successfully fit not only the / = S-wave, but also reproduce the repulsive 
1 — 2 S-wave. Their pole positions are given in Table 1. 

Schechtcr's group has also contributed a series of papers on all of cr, k, 
/o(980) and ao(980) [17-21]. This work examines possible mixing between 
2-quark and 4-quark states. 

Zheng and collaborators at Peking University have developed new types 
of dispersion relations and have applied them to analysis of data on the cr 
and n [22-25]. In particular, they clarify how unphysical sub-threshold poles 
can originate from the 1/v^ dependence of phase space factors; they also 
show that a pole is possible without the phase shift reaching 90°, as may well 
happen in Kti scattering. Thirdly, in Ref. 25, they repeat the analysis of 
Colangelo et al. with different dispersion relations fitting the left-hand cut; 
they confirm the necessity for the a pole. 

6 K-matrix fits * 

Focus data on Svr exhibit a similar low mass peak to E791 data [26]. If 

it is fitted with a Breit-Wigner with T oc p^Tr, there is a pole at M = (455 ± 
36) — i(190±36) MeV [27]. However, Focus show that their data can also be 
fitted with a K-matrix parametrisation of Anisovich and Sarantsev [28] which 
does not include a a pole. This provides an escape route for those who do 
not wish to believe in the a pole. The Babar collaboration has likewise found 
that their data can be fitted by the same K-matrix parametrisation [29] . This 
point requires some detailed and (unfortunately) critical comments. 

There are several problems. The K-matrix parametrisation includes low 
mass poles at ~ 600 and 1200 MeV and at s = — 3 GeV^. Focus report their 
parametrisation fully, so it is straightforward to follow details of the fit. 

The first point is that the low mass region is being fitted by 3 poles instead 
of 1. This introduces large flexibility into the fit. 

The second point is that the pole below threshold plays a major role, 
but is questionable. If it is a zero of -D(s), it is a very deeply bound state 
which is more questionable than the a. A more likely interpretation is a 
pole of N(ys). However, the N function is different for tttt elastic scattering 
and production reactions. So there is no reason to think that a pole fitted 
to elastic scattering will be appropriate to production. Furthermore, the N 
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function is normally interpreted as providing driving forces which generate 
the amplitudes on the right-hand cut, e.g. resonances. Including them as 
specific poles in searching for resonances looks like double-counting. 

There is a third problem. In the K-matrix approach, tttt elastic scattering 
below the KK threshold is fitted not only by Kn (i.e. tttt tttt) but also 
by the analytic continuation of K12 below the KK threshold. In the case of 
D decays, this corresponds to D ^ {KK)n, followed by KK rescattering to 
TTTT. This component should be constrained to reproduce data on — > KKtt 
above threshold; however, that has not been done. 

In K-matrix fits to tttt elastic scattering below the KK threshold, I know 
from personal experience that there is large cross-talk between Ku and K12. 
K12 increases below threshold as \k\/^/s where \k\ is the magnitude of the 
virtual momentum in the KK channel. The increasing K12 as s — > looks 
not unlike a a pole. The extrapolation below KK threshold is reliable only 
close to threshold. In the Flatte formula for /o(980), as an example, it 
plays an important role in contributing to the real part of the amplitude 
near resonance. However, far from threshold the extrapolation is hazardous. 
This is a general problem which has been seen elsewhere. For example, in 
calculations of charmonium levels, open charm states with large decay widths 
perturb energy levels by major amounts. 

My view is that there is a stabilising factor which needs to be brought into 
play. It is common experience that form factors play a role above threshold 
with a radius of order 0.6-0.8 fm. If such a form factor is needed above 
threshold, the analytic continuation below threshold breaks down because 
it depends critically on the imaginary part of the amplitude extending to 
infinity. With a form factor, the continuation below threshold can instead 
be done using a dispersion integral: 



I have carried out calculations along these lines. The form factor gives Im K12 
a localised peak close to threshold and the result is that Re K12 drops fairly 
rapidly below threshold. In practice, it falls below threshold with a form 
similar to the fall-off above threshold, e.g. exp(— aj/cp). This behaviour is 
totally different from the analytic continuation without a form factor, and 
can easily lead to an order of magnitude reduction in Re K12 by the time 
one reaches the position of the a pole. Focus show in their Figs. 7 and 8, the 
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magnitudes they fit to sub-threshold KK^ An, rjrj and 7777' amphtudes; all are 
very large compared with fitted values of these amplitudes above threshold, 
typically by one to two orders of magnitude. It seems physically unreasonable 
that channels like r^r] and ryr;', which are quite small in the physical region, 
should make large sub-threshold contributions. 

Putting this point in a different way, why stop at the 47r channel above 
threshold? Why not include an infinity of open channels at high mass: Gtt, 
Stt, ... IOOtt . . . This leads to the possibihty of fitting all mesons resonances 
as interferences between an infinity of distant high-mass singularities. This 
contradicts the conventional view that low mass peaks may be attributed to 
nearby singularities unless there is a specific convincing alternative. 

My impression is that the K-matrix is dependable above threshold, where 
it is fitted directly to data on the inelastic channels. But below threshold, it 
introduces excessive fiexibility. What is urgently needed is to constrain it to 
fit the left-hand cut as well, following the work of Colangelo et al. 

Conversely, a criticism of published fits based on the T-matrix is that they 
ignore sub-threshold contributions from 47r, KK and r/r/. A more rigorous 
expression for the Breit-Wigner denominator of eqn. (1) is 



In work not yet published, I have used this full form. It takes account of 
the opening of the inelastic channels and also the effect of the dispersive 
contribution in eqn. (3) to Re Ku. The detailed parameters of the fit 
change, but effects on the position of the a pole are within the presently 
quoted errors. The essential reason for this is that the BES data define the 
position and width of the peak unambiguously. What happens is that the 
additional term m{s) in the denominator is accomodated by corresponding 
changes in the already fiexible form of N[s). 

7 The K pole 

E791 data on — > (ir~7r+)7r+ provided the first evidence for the k, pole 
from production data [30]. Their parametrisation gives a pole at M = (721 ± 



D{s) 



- s - m{s) - iNtot{s) 



(6) 
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Figure 6: BES II data for J/\l/ K^K^tt^ti^ . (a) The scatter plot 
M(K^Ti^) V. M{K~7i^). Projections on to (b) K'^tt^ mass, (c) and (d) 
vrvr and ttttK mass, (e) The Dalitz plot for events where one K^ir^ pair has 
mass 892 ± 100 MeV. (f) Mass projection of the second K^tc^ pair for the 
same selection as (e). 
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Figure 7: (a) The mass projection of Fig. 6(f) divided by Kti phase 
space, in bins of 20 MeV. 

61) — i(292±131) MeV. A combined fit to these and other data will be shown 
below. 

Next, BES II data on J/^f ^ K^ti~K~tt^ provide further evidence for 
the K in the channel J/^ iir*(890)fi; [31,32]. The scatter plot and mass 
projections are shown in Fig. 6; histograms show the fit. There are clear 
peaks due to K*{%9Q), ^^0,2(1430), p(770), p(1270), 7^1(1270 + 1400) and 
7^1(1770). If one selects K^-k"^ pairs within 50 MeV of 892 MeV, the mass 
projection of the other i^^Tr^ pair is shown in panel (f). A broad i^vr S-wave 
component is visible under the ir*(890). Fig. 7 shows the effect of factoring 
out Xtt phase space in the 4-body system. There is a broad peak below 750 
MeV, which is fitted as the n signal. 

Data over the whole of phase space are fitted. Full details are given in 
Ref. [32]. The fit is made simultaneously to BES data and LASS data [33] 
on the Kt^ I — 1/2 S-wave. Eqn. (1) is used with 



with sa — rn\ — 0.5m^. [Fits of similar quality may be obtained by re- 
placing exp(— ay^) with exp(— a's) or l/(s — sq) with marginally poorer log 
likelihood]. If the factor (s — s a) is omitted from eqn. (8), the poor fit is 
shown in Fig. 8(a). The k mass projection is shown by the full histogram of 



N{s) = M{s - Sa) ex.p{-ay/s)pKn{s) 



(8) 
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Figure 8: (a) The poor fit when K*{890)k is removed; (b) individual con- 
tributions from K (full histogram) and ii'o(1430) (dashed); (c) the coher- 
ent sum oi K + ii'o(1430) (dashed histogram) and the coherent sum k + 
Ko{U30) + /s:i(1270) + i^'i(1400) (full). 



Fig. 8(b) and the -^"0(1430) mass projection as the dashed histogram. There 
is destructive interference between k and /('o(1430); their coherent sum is 
shown by the dashed histogram in Fig. 8(c). Sensitivity to this interference 
is one reason for fitting the LASS data simultaneously; it is absent in those 
data, where the phases of the two components add. A second reason is to 
examine whether eqn. (8) fits both sets of data successfully. 

Since Ki decays populate the low mass Ktt range, it is essential to demon- 
strate that Ki decays do not reproduce the k peak. Omitting the k leads to 
a fit worse by > 1000 in log likelihood. 

An essential feature of the analysis is the separation of the 1430 MeV 
Ktt peak between ^^0(1430) and i^2(1430). This separation can be made 
cleanly only by analysing angular correlations between i^'*(1430) and the 
accompanying Kj{lA30) ('entanglement'). The result is that the peak is 
(75 ± 3)%KoiU30) and 25%ir2(1430). 

The fi pole optimises at M = (760 ± 20 {stat) ±AO{syst))-i (420 ± 45 ± 60) 



14 



MeV. The Ko(1430) is fitted with a Flatte formula including coupling to Krj'-^ 
it requires g'^{Kr]')/ g'^{Kr]) = 1.0 ± 0.3. A detail is that the Adler zero is 
also included into the width of the Kq{1A30) so that there is an Adler zero 
in the full Ktt S-wave; this zero improves the fit noticeably. 

8 The tv phase 

The phase variation of the k with mass is well determined in two ways which 
agree. Firstly, there is a large interference between channels K*{980)k and 
-ft'fri(1270 + 1400). Secondly, there is a large interference between k and 
-^"0(1430), which both contribute to the Ktt S-wave. 

A bin- by-bin fit has been made where the k signal is fitted in magnitude 
and phase in 10 individual bins 100 MeV wide. Results are shown in Fig. 9 
and are in good agreement with the global fit (full curves). [The data points 
have been adjusted for the phase (f)o of the isobar model, so that the k phase 
is zero at threshold]. The agreement between the bin- by-bin fit and the full 
curve shows that the phase variation can be fitted by the same D{s) in both 
cases, i.e. purely by a resonance. 
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Figure 9: Points show (a) the phase (b) the magnitude of the k amplitude, 
determined bin-by-bin; curves shows the global fit. 

There is a point of general interest concerning the k phase. If there is 
a pole near threshold, why does the observed phase not go through 90° at 
low mass? This is a subtle point. The clue is that the k, pole lies in the 
complex s-plane almost below the KK threshold. The factor {s — sa) in 
eqn. (8) develops an imaginary part as one goes off the real s-axis. Also 
the phase space factor p{s) becomes complex. The combined effect of these 
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two effects is to produce a phase rotation of nearly 90° as one moves off the 
real s-axis to the pole. At the pole, the phase does indeed go through 90°. 
But on the real s-axis, the constraint of unitarity holds the phase at zero 
at the Kn threshold. There is a retardation of the k. phase along the real 
s- axis of roughly 90° with respect to the pole, although there is also some 
s-dependcncc. This was pointed out by Oiler [3]; a corresponding smaller 
effect occures likewise for the a pole. On the real s-axis, one is really seeing 
the upper side of the k and a poles but with retarded phases caused by the 
curious effects of analyticity near threshold. This is an unfamiliar situation 
and has caused widespread confusion. 

A criticism has been made that n and iro(1430) cannot be separated in 
the Kn S-wave. If so, phases could be adjusted freely in the bin-by-bin fit, 
but the fit resists that. Secondly, the criticism is logically incorrect. The ar- 
gument is that magnitudes and phases of both k and ii'o(1430) can be fitted 
freely as a function of s. If that were true, it would always be impossible 
to separate any two resonances with the same in a single channel. How- 
ever, analyticity requires that real and imaginary parts of the amplitude are 
related and cannot be fitted independently. In the conventional approach, 
each resonance is fitted by a Breit-Wigner amplitude where magnitude and 
phase are parametrised by a single function of s. There is indeed flexibility 
in N{s) for the n, but only limited flexibility; to a first approximation, the 
real part of the amplitude is close to the gradient of the imaginary part, as 
in the simple Breit-Wigner formula. The upshot is that the observed inten- 
sity in the Kir S-wave as a function of s is determined by the intensities of 
the individual k and 7^'o(1430) and the real part of the interference between 
them. These three terms have distinctively different s-dependence and can- 
not be confused, except within statistics of each component. Furthermore, 
interferences of both k and Ko{U30) with KKi(1270 + 1400) provide an in- 
dependent determination of phases. The agreement with the phase in LASS 
data also seems hardly fortuitous. 

9 The BES fit to the k 

An independent analysis of exactly the same data has been reported by the 
BES group [31]. I do not have the figures, so it is necessary to refer to 
Ref. [31]. This analysis parametrises the k in the same general form as 
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the Ishida model [34]. The Kn S-wave amphtude is fitted as the sum of 
a conventional Breit-Wigner amplitude with F oc pxir plus an incoherent 
background which is adjusted to fit the remaining intensity as a function of 
mass. The Breit-Wigner amplitude is fitted with M = 878 ± 23ttt MeV 
with r at this mass of 499 ± 45I72 MeV. The corresponding pole position is 
M = (841 ± SOtjl) - z(309 ± 451^^ MeV. 

There are some important technical differences between the two analyses. 
In my analysis, one of the phase determinations comes from interference be- 
tween K and iiro(1430). Ref. [31] states that the BES analysis omits this inter- 
ference. Secondly, because of the interference between ii"o(1430) and k, it is 
important to separate the channels ir*(890)i^2(1430) and K*{890)Ko{U30). 
This separation requires analysis of the full angular correlations between de- 
cays of ir*(890) and Kj{U30). The BES analysis does not consider K*{890) 
decays [31] and therefore separates -^"0(1430) and -ft'2(1430) only from the 
line-shape. Thirdly, the BES analysis does not include interference between 
K and Ki —>■ pK. For these reasons, the phase of the k is not determined in 
that analysis. 




Figure 10: The Dalitz plot for E791 data for D+ {K-n+)n+. 
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10 Re-analysis of E791 data 



The Dalitz plot for E791 data on D+ {K-n+)n+ is shown in Fig. 10. 
Along the 7^'*(890) bands, there is an obvious asymmetry due to interference 
with the Kn S-wave. A new fit has been reported recently where the magni- 
tude and phase of the Kn S-wave anplitude are fitted separately in 37 mass 
bins [35] . Results for the magnitude and phase of the Ktt S-wave amplitude 
will be shown below in Figs. 11 and 12(f). 

I have carried out a combined fit to these new E791 data, together with 
LASS data and BES, with the objective of separating k, and 7i'o(1430) signals. 
The BES data define well the i^ro(1430) peak, which is much more conspicuous 
there than in either LASS or E791 data. To a first approximation, the BES 
data do not determine the k parameters strongly; they are determined mainly 
by the LASS and E791 data. 

In the E791 fit (and also Ref. [30]), the amplitude for production includes 
a form factor F = exp(— ag^), where q is the momentum of the k in the Ds 
rest frame and a = 2.08 GeV~^. I have varied the exponent a and results are 
shown in Fig. 11. Panel (a) uses a — and achieves the best fit. In (b)-(d) 
a is increased in equal steps to the E791 value in (d). There is an obvious 
preference for a close to 0. In fact a optimises at within experimental 
error for both magnitude and phase. This corresponds to a point-like decay 
Dg — > KTT with an RMS radius < 0.38 fm with 95% confidence. 

Fits to the BES II data and LASS data are shown in Fig. 12. In BES 
data, the i^o(1430) is a large signal with well defined centroid and width. 
The fit to this peak is shown in Fig. 12(e) using 25% K*{890)K2{U30) and 
75% K*{890)Ko{U30), as determined from the fit to BES data alone. Figs. 
12(c) and (d) show that the fit to the magnitude and phase of the k in BES 
data is acceptable. Figs. 12(a) and (b) show the fit to LASS data. There 
is a small systematic discrepancy around 1.2 GeV which resists a variety of 
parametrisations. It could be associated with the onset of KKmr inelasticity 
which is presently unknown. 

From the combined fit, the k pole position is M = (750;55) - i(342 ± 60) 
MeV. This compares with a fit to LASS data alone M = (722±60) -i(386± 
50) MeV [36]. This interpretation of E791 data brings results into good 
agreement with both LASS and BES II data. 
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Ktt mass (GeV) Ktt mass (GeV) 

Figure 11: Fits to the magnitude of the k amphtude in E791 data from Ref. 
[35], for four values of a in the form factor; a is in units of GeV^^. 

11 /o(980) 

BES II data on J/^' — > (j)7T~^n~ and (j)K^K~ both contain prominent /o(980) 
signals [37]. The data are fitted with the Flatte formula 

/ = 1/[M^ - s - i{glp^^ + qIpkk)]. (9) 

Fitted parameters are M = 965 ± S{stat) ± Q{syst) MeV, gl = 165 ± 10 ± 15 
MeV, gl/gl = 4.21 ± 0.25 ± 0.21. There is a second sheet pole at (998 ± 4) - 
i(17 d= 4) MeV, very close to the KK threshold, and a distant third sheet 
pole at (851 ± 28) — 2(418 ± 72) MeV. The dominance of the narrow second 
sheet pole is used by Baru et al. [38] to argue that the /o(980) is mostly a 
KK bound state resembling the deuteron. 

The best current determination of the parameters of ao(980) is from Ref. 
[39]. The second sheet pole is at M = (1032 ± 5) - i(85 ± 9) MeV, and the 
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Kn mass (GeV) Kn mass (GeV) 



Figure 12: Fits to (a) and (b) the phase and magnitude of LASS amplitudes 
for elastic scatering, (c) and (d) the phase and magnitude in BES data, (e) 
the 1430 MeV peak in the Kn mass projection of BES data, (f) the fit to 
E791 phases with a = 0. 



third sheet pole at M = (981 ±8) -i(304±26) MeV. This behaviour is much 
more like a conventional resonance than /o(980). 

It is often argued that they are anomalously narrow states and therefore 
decay weakly. This is quite wrong. The couphng of /o(980) to KK corre- 
sponds to a width of 695 MeV when pxK ^ 1- The resonance appears as a 
narrow cusp at the KK threshold because the KK channel opens extremely 
rapidly; the increase in the total width pulls the intensity in the tttt channel 
down sharply at threshold, making the resonance appear narrow. 

There is another interesting way of viewing this effect. The amphtude for 



20 



KK — > TTTT is given by 

1 '^g-Kg-K\lKkK/s 



(10) 



where /c^r and /ck are centre of mass momenta. The cross section for this 
process goes as 1//ck near threshold; this is the famihar 1/v law for decay to 
open channels near threshold. The imaginary part of the amplitude follows 
the optical theorem 

Im f = kKCTtot (11) 

and therefore has a step at threshold. Prom the dispersion relation for the 
amplitude 

^ l.W(^ (12) 

the effect of the step is to generate a narrow peak in the real part at threshold. 
This is illustrated in Fig. 13 for the scattering length observed for /o(980). 

It is a general feature of the opening of any threshold. It provides a mecha- 
nism by which a resonance can be attracted to a threshold by the additional 
attraction of the threshold. It seems likely that this is the explanation of 
how /o(980) and ao(980) are captured by the KK threshold. 




Figure 13: Real and imaginary part of the amplitude for KK — > tttt. 
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12 What are a, /o(980) and ao(980)? 



12.1 (7 models and Chiral Symmetry breaking 

The saga begins in the early 1960's. At this time, it was a puzzle why 
the lightest mesons were the tt and the K, whereas ground-states of even- 
even nuclei were 0+. This led Schwinger to postulate the existence of a cr 
meson. A second puzzle was the Goldberger-Treiman relation and the fact 
that the axial weak current is nearly conserved, like the vector current. In an 
effort to save the Goldberger-Treiman relation from renormalisation effects, 
Gell-Mann and Levy invented Chiral Symmetry [40] and the linear a model, 
placing the pion and the cr on an equal basis. Nambu and Jona-Lasinio 
invented the non-linear a model [41]. Today it is still an open question 
whether these models have any connection with the observed a pole. 

When QCD emerged, it was clear that Chiral Symmetry is broken by the 
appearance of mass terms of the form mqq in the Lagrangian. Today it is 
almost universally accepted that Chiral Symmetry is spontaneously broken. 
Spontaneous breaking of a symmetry leads to the existence of a massless 
Goldstone boson. Many examples are known in condensed matter physics. 
A classic example is ferromagnetism where, below a critical point, spins align 
spontaneously; there is an associated massless magnon responsible for spin- 
waves. In a crystal, the regular spacing of atoms spontaneously creates order; 
the phonon is the associated massless particle. 

It is widely believed that the pion is almost massless because of Chiral 
Symmetry breaking; this leads to the Adler zero, which has played a promi- 
nent role in fitting the data on cr and ac. The small mass of the pion derives 
from the small masses of u and d quarks. Gell-Mann, Oakes and Renner 
showed [42] that 



Here, f^^ is the pion decay constant. It is intriguing that the mass of the pion 
is related to the density of 0"*" fluctuations. It is not presently clear whether 
these fluctuations have any relation to cr, etc. 

A conflned quark will surround itself with a cloud of gluons and sea- 
quarks. Two sets of authors [43,44] suggest that the resulting constituent 
quarks will acquire a mass of roughly one-third the mass of the nucleon. A 
cr made of two such constituent quarks would have roughly the right mass. 




(13) 
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However a problem with this scheme is that one would then expect the a 
to have a brother with / = 1 at a very similar mass, whereas the ao(980) is 
~ 400 MeV heavier. 

12.2 Jaffe's model 

This problem led Jaffa to propose that a and its relatives are q^q^ states 
[45] . His suggestion is that there is a pairing interaction forming diquarks in 
the flavour 3 configuration: ud, ds and us. Then 3 and 3 make a colourless 
nonet. The a is the / = member uddu, the /t"*" is usdd, aj] is ss{uu — 
and /o(980) is ss{uu + dd)/\/2. This scheme neatly explains why oq and 
/o are degenerate in mass and heavier than the a by twice the mass of the 
s-quark. It also neatly fits in with the intermediate mass of the k. Maiani et 
al. advocate the same scheme [46]. 

An interesting experimental question is whether there is a further ssss 
state. This is foreign to Jaffe's scheme. GAMS have reported tentative 
evidence for a narrow state in rjr]' at 1914 MeV, almost exactly the rj'r)' 
threshold [47]. Such a state would decay easily to rjrj' and would be more 
prominent there than in 7]'r]', just as /o(980) is more prominent in J/^ — > 
07r"'"7r~ than in J/^ — > (f)K~^K~: more phase space. 

There is support for Jaffe's scheme from recent Lattice QCD calculations 
of Okiharu et al. [48]. They find configurations of Fig. 14(a) at large radii 
and those of Fig. 14(b) at small r. Lattice QCD calculations favour the 
multi- Y-shaped flux-tube configuration for the connected 4-quark system. 
One can rationalise this scenario with the argument that quarks at large r 
are non-perturbative and acquire mass from dressing. Another view of this 
is that it costs energy to establish a flux tube between quarks separated 
radially. The massive g^g^ configuration can decay by fission to two hghter 
pions at small r. 

There is experimental support for the notion that decays to irn are at 
short range. The a amplitude is known v. momentum k up to nearly 2 
GeV/c. Its Fourier transform then determines the radial dependence of the 
matrix element; the RMS radius of the matrix element is 0.4 fm [36]. 

There is one important problem with Jaffe's scheme in its simplest version. 
The ratio r = g^{fo{980) K K) / g\ao{980) KK) = 2.7 ± 0.5 disagrees 
with the ratio 1 predicted from Jaffe's model. [One cannot escape from this 
problem if /o(980) and ao(980) are 2-quark states. Detailed arithmetic on 
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Figure 14: (a) a connected tetraquark configuration, (b) a 'two-meson' con- 
figuration. 

branching ratios, making use of mixing between /o(980) and cr, predicts a 
ratio sliglitly below 1.] 

Tornqvist [49] points out that both ao(980) and /o(980) must contain 
a KK component in their wave functions. His eqn. (15) includes mesonic 
channels AiBj, in loop diagrams: 



Here Re n(s) = g^j^y4m^/s— 1 for s < 4m|;. [Tornqvist's equation is 

written in terms of qq, but could equally well be reformulated in terms of 
4-quark states]. His formula is easily evaluated to find the KK components 
in ao(980) and /o(980) as functions of s. At the KK threshold, the binding 
energy — > and the KK wave function extends to infinity, so the KK fraction 
— > 1. Results are shown in Fig. 15 by the dotted curves. This figure also 
shows hne-shapes as the full curves. The mean KK fraction integrated over 
the hne-shape is > 60% for /o(980) and half this for ao(980). 

Suppose decays to nn and ttt] occur at small r, following the Lattice QCD 
picture. [There will also be some 4g KK decay at small r]. When such 
decays occur, the passive kaonic cloud at large r is left 'in the air' (adiabatic 
approximation), and contributes to 'fall-apart' decays. The intensity of these 
fall-apart decays is proportional to the KK intensity in the wave function, 
and reduces the discrepancy with Jaffe's model by a factor 2. However, this 
argument docs not account for /o(980) lying closer to the KK threshold than 
ao(980); that must be taken from experiment. 



> +E.[-(^//^/.s)Re ni(s)]^-/^|A,Di > 

1 - j:i{d/ds)Re ni(s) 



(14) 




24 



: a)fo(980) / 




_KK "Taction j 


Ine-shope 

\ / 


^ R(fm)/10 ^ — - 


0.9 


1 1 




M(7T7t) (GeV) 



: b) ao(980) / 

/ 1 


\ 


' KK f Taction ./ ■■; 

- \ ' 


\ 1 i r e — s h □ p e 


l-^ R(fm)/10 


0.9 


1 1 

M(7T7i) (GeV) 



Figure 15: Line shapes of /o(980) and ao(980), full curves; the KK fraction 
in the wave function, dotted curves; O.lxRMS radius, dashed curves. 



12.3 CDD poles 

There is another important distinction between a (and its family) and the 
accepted qq states. In tttt elastic scattering, p and /2 exchanges in u and 
t channels generate attraction. If one takes the X-matrix in the s-channel 
from these Born terms, the unitarised amplitude K/{1 — ipK) reproduces 
the observed tttt S-wave quite well up to 1 GeV and beyond [50,51]. Similar 
exchanges are roughly sufficient to generate /o(980) and ao(980). 

In this respect, a and its relatives behave completely differently to regular 
qq states like p, ui, i^*(890) and (p. These cannot be derived from t and u- 
channel exchanges. It was this fact which led to the quark model. The 
commonly accepted qq resonances appear as CDD poles [52] ; these are poles 
which identically satisfy dispersion relations with no apparent connection 
with driving terms in N{s) on the left-hand cut. 

For this reason, some people view a, k, /o(980) and ao(980) as molecular 
states made of mesons. There is, however, an alternative view. In nucleon- 
nucleon physics, the observed partial waves can be explained in terms of 
meson exchanges. It is conventional to view these exchanges as generating 
a potential V{r). In meson physics, it is equally possible to take the view 
that Tin elastic scattering is telling us about the confining potential itself. 
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Regular qq states are confined within this potential; their leakage out of the 
potential dictates their decay widths. These widtlis must come out so that 
they satisfy analyticity over left and right-hand cuts. 

Any connection of the a, k, /o(980) and ao(980) poles to this potential is 
presently speculative. The formation of nuclei is clearly a phase transition. 
Chiral Symmetry breaking is also clearly a phase transition; Lattice QCD 
calculations tell us that the phase boundary is at about 160 MeV. An impor- 
tant question if whether Confinement is the same phase transition as Chiral 
Symmetry breaking. Present indications from Lattice QCD calculations are 
that it has a similar transition temperature. 

12.4 The scheme of Van Beveren and Rupp 

Another interesting, related scheme is that of Van Beveren and Rupp [53-55] . 
They set out to model the spectrum of all mesons from the lightest to char- 
monium and bottomonium. This is done with either a harmonic oscillator 
potential [53], which can be handled algebraically, or any kind of confine- 
ment spectrum in a more general approach [54,55]. To allow for decays, they 
couple qq states to outgoing mesons through a transition potential. Ideally, 
this transition potential has a ^Pq dependence on radius r [53], but in most 
calculations it is approximated by a 5 function at r = 0.6-0.7 fm or smaller 
(for heavy q-qbar systems). Their equations include a relativistic reduced 
mass for the two outgoing mesons. This reduced mass generates a zero in 
the amplitude very similar to the Adler zero. 

In this scheme, /o(1370), ao(1450), A'o(1430), etc. are regular though 
unitarised qq states; a, k, /o(980) and ao(980) appear as 'extra' states created 
by coupling of qq to mesonic decays. The a and n were indeed predicted 
rather well in 1986 [53]. It is instructive to vary the coupling constant A 
between confined states and mesons. For very weak coupling, a and k appear 
as continuum states with infinitely large widths, i.e. at very large —Im y^. 
Fig. 16 shows the analogous Ds{2317) situation. As the coupling between 
confined qq states and outgoing mesons is increased, the continuum state 
moves towards the real axis; simultaneously, the cs state acquires a width 
and moves off the real s axis. In this way, they account for the Ds{2317) 
like (7 and k, i.e. as an extra, "molecular" state, created by coupling of the 
DK continuum to regular cs 0+ state; the latter are pushed to a higher mass 
than that of the usual funnel potential [55]. 
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Figure 16: Trajectories as a function of coupling constant A to decays. 
12.5 Extrapolation to 3 

Theorists have tried varying the numbers of colours using Unitariscd Chi- 
ral Perturbation Theory [56,57]. They find that as Nc increases, a and n 
fade away into the continuum. [Ref. [58] agrees that they change by large 
amounts, but reaches different conclusions about the precise change]. This 
is in contrast to the p which survives largely unchanged as Nc increases. 
That is further evidence that the light scalars have a different character to 
conventional qq states. 
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